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' Abstract 

I The general relativistic Lense-Thirring precessions of the periheHa 

of the inner planets of the Solar System are < 10~^ arcseconds per cen- 
^ ' tury. Recent improvements in planetary orbit determination may yield 

the first observational evidence of such a tiny effect. Indeed, corrections 
to the known perihelion rates of -0.0036±0.0050, -0.0002±0.0004 and 
1^ . 0.0001±0.0005 arcseconds per century were recently estimated by E.V. 

' Pitjeva for Mercury, the Earth and Mars, respectively, on the basis of 

. the EPM2004 ephemerides and a set of more than 317,000 observations 

of various kinds. The predicted relativistic Lense-Thirring precessions 
O . for these planets are —0.0020, —0.0001 and —3 x 10~^ arcseconds per 

century, respectively and are compatible with the determined perihelia 
iH^. corrections. The relativistic predictions fit better than the zero-effect 

hypothesis, especially if a suitable linear combination of the perihe- 
lia of Mercury and the Earth, which a priori cancels out any possible 
bias due to the solar quadrupole mass moment, is considered. How- 
ever, the experimental errors are still large. Also the latest data for 
' Mercury processed independently by Fienga et al. with the INPOP 

ephemerides yield preliminary insights about the existence of the solar 
Lense-Thirring effect. The data from the forthcoming planetary mis- 
sion BepiColombo will improve our knowledge of the orbital motion 
of this planet and, consequently, the precision of the measurement of 
the Lense-Thirring effect. As a by-product of the present analysis, it 
is also possible to constrain the strength of a Yukawa-like fifth force to 
a 10^^^ — 10^^^ level at scales of about one Astronomical Unit (10^"'^ 
m). 

Keywords: experimental tests of gravity, Lense-Thirring effect, orbital 
motions, planets, Solar System 
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1 Introduction 

A satisfactorily empirical corroboration of a fundamental theory requires 
that as many independent experiments as possible are conducted by different 
scientists in different laboratories. Now, the general relativistic gravitomag- 
netic Lense-Thirring (LT) effect is difficult to test, also in the weak-field and 
slow-motion approximation, valid in our Solar System, both because such a 
relativistic effect is very small and the competing classical signals are often 
quite larger. Until now, a 6% LT test has recently been conducted in the 
gravitational field of Mars by using the data of the Mars Global Surveyor 
(MGS) spacecraft; other tests accurate to abouil] 10% have been performed 
by a different team in the gravitational field of the Earth by analyzing the 
data of the LAGEOS and LAGEOS II artificial satellites (see Section [1.2p . 
We think it is worthwhile to further extend these efforts trying to use differ- 
ent laboratories, i.e. other gravitational fields, even if the outcomes of such 
tests should be less accurate than those conducted so far. To this aim, in 
this paper we investigate the LT effect induced by the Sun on the orbital 
motion of the inner planets of the Solar System in the context of the latest 
results in the planetary ephemerides field. 

1.1 The Lense-Thirring effect 

According to Einstein, the action of the gravitational potential [/ of a given 
distribution of mass-energy is described by the coefficients gfj,^, /i, = 
0,1,2,3, of the space-time metric tensor. They are determined, in prin- 
ciple, by solving the fully non-linear field equations of the Einsteinian Gen- 
eral Theory of Relativity (GTR) for the considered mass-energy content. 
These equations can be linearized in the weak-field {JJ /(? << 1, where c is 
the speed of light in vacuum) and slow-motion {v/c « 1) approximation 
(Mashhoon 2001; Ruggiero and Tartaglia 2002), valid throughout the Solar 
System, and look like the equations of the linear Maxwellian electromag- 
netism. Among other things, a noncentral, Lorentz-like force 




(1) 



Other estimates point towards a 15 — 20% error. 
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acts on a moving test particle of mass m. It is induced by the post- 
Newtonian component Bg of the gravitational field in which the particle 
moves with velocity v. Bg is related to the mass currents of the mass- 
energy distribution of the source and comes from the off-diagonal compo- 
nents goi,i = 1,2,3 of the metric tensor. Thanks to such an analogy, the 
ensemble of the gravitational effects induced by mass displacements is also 
named gravitomagnetism. For a central rotating body of mass M and proper 
angular momentum L the gravitomagnetic field is 

G[3r(r ■ L) -r^L] 

cr^ 

One of the consequences of eq. ([T]) and eq. ([2]) is a gravitational spin- 
orbit coupling. Indeed, if we consider the orbital motion of a particle in the 
gravitational field of a central spinning mass, it turns out that the longitude 
of the ascending node O and the argument of pericentre to of the orbit 
of the test particle are affected by tiny secular advances J^lt, i^lt (Lense 
and Thirring 1918, Barker and O'Connell 1974, Cugusi and Proverbio 1978, 
Soffel 1989, Ashby and Alhson 1993, lorio 2001) 

• 2GL . 6GLcosi 

"LT = 3", f^LT = J, (3) 

where a, e and i are the semimajor axis, the eccentricity and the inclina- 
tion, respectively, of the orbit and G is the Newtonian gravitational con- 
stant. Note that in their original paper Lense and Thirring (1918) used the 
longitude of pericentre vj = Q + uj. 

The gravitomagnetic force may have strong consequences in many astro- 
physical and astronomical scenarios involving, e.g., accreting disks around 
black holes (Thorne et al. 1986; Stella et al. 2003), gravitational lensing and 
time delay (Sereno 2003; 2005a; 2005b). Unfortunately, in these contexts 
the knowledge of the various competing effects is rather poor and makes 
very difficult to reliably extract the genuine gravitomagnetic signal from the 
noisy background. E.g., attempts to measure the LT effect around black 
holes are often confounded by the complexities of the dynamics of the hot 
gas in their accretion disks. On the contrary, in the solar and terrestrial 
space environments the LT effect is weaker but the various sources of sys- 
tematic errors are relatively well known and we have the possibility of using 
various artificial and natural orbiters both to improve our knowledge of such 
biases and to design suitable observables circumventing these problems, at 
least to a certain extent. 
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1.2 The performed and ongoing tests 

Up to now, all the performed and ongoing tests of gravitomagnetism were 
implemented in the weak-field and slow-motion scenarios of the Earth and 
Mars gravitational fields. 

As far as the Earth is concerned, in April 2004 the GP-B spacecraft 
(Everitt 1974; Fitch et al. 1995; Everitt et al. 2001) was launched. Its aim 
is the measurement of another gravitomagnetic effect, i.e. the precession 
of the spins (Pugh 1959; Schiff 1960) of four superconducting gyroscopes 
carried onboard. The level of accuracy obtained so far is about 256 — 128% 
(Muhlfelder et al. 2007), with the hope of reaching ~ 13% in December 
2007. 

In regard to the LT effect on the orbit of a test particle, the idea of 
using the LAGEOS satellite and, more generally, the Satellite Laser Ranging 
(SLR) technique to measure it in the terrestrial gravitational field with the 
existing artificial satellites was put forth for the first time by Cugusi and 
Proverbio (1978). Attempts to practically implement such a strategy began 
in 1996 with the LAGEOS and LAGEOS II satellites (Ciufolini et al. 1996). 
The latest test was performed by Ciufolini and Pavlis (2004). They analyzed 
the data of LAGEOS and LAGEOS II by using an observable independently 
proposed by Pavlis (2002), Ries et al. (2003a, 2003b) and lorio and Morea 
(2004). The error claimed by Ciufolini and Pavlis (2004) is 5-10% at 1-3 
sigma, respectively. The assessment of the total accuracy of such a test 
raised a debate (lorio 2005a; 2005b; 2006a; 2007; Ciufolini and Pavhs 2005; 
Lucchesi 2005). 

Recently, a 6% LT test on the orbit of the Mars Global Surveyor (MGS) 
spacecraft in the gravitational field of Mars has been reported (lorio 2006b); 
indeed, the predictions of general relativity are able to accommodate, on 
average, about 94% of the measured residuals in the out-of-plane part of the 
MGS orbit over 5 years. 

Finally, it must be noted that, according to Nordtvedt (2003), the multi- 
decade analysis of the Moon'orbit by means of the Lunar Laser Ranging 
(LLR) technique yields a comprehensive test of the various parts of order 
0{c~'^) of the post-Newtonian equation of motion. The existence of the 
gravitomagnetic interaction as predicted by GTR would, then, be inferred 
from the high accuracy of the lunar orbital reconstruction. A 0.1% test 
was recently reported (Murphy et al. 2007): a critical discussion of the real 
sensitivity of LLR to gravitomagnetism can be found in (Kopeikin 2007). 

^Se e on the W EB StanfordNews 4/14/07 downloadable at 
|http://einstein.stanford.edu/[ 
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Table 1: Gravitomagnetic secular precessions of the longitudes of perihelion 
w of Mercury, Venus, Earth and Mars in " cy~^. The value (190.0 it 1.5) x 
10^^ kg s~^ (Pijpers 1998; 2003) has been adopted for the solar proper 
angular momentum Lq. 

Mercury Venus Earth Mars 

-0.0020 -0.0003 -0.0001 -3 x 10"^ 



Also the radial motion of the LAGEOS satellite would yield another indirect 
confirmation of the existence of the gravitomagnetic interaction (Nordtvedt 
1988). 

2 The solar gravitomagnetic field 

The action of the solar gravitomagnetic field on the Mercury's longitude of 
perihelion was calculated for the first time by de Sitter (1916) who, by as- 
suming an homogenous and uniformly rotating Sun, found a secular advance 
of —0.01 arcseconds per century ( " cy~^ in the following). This value is also 
quoted by Soffel (1989). Cugusi and Proverbio (1978) yield -0.02 " cy"^ for 
the argument of perihelion of Mercury. Instead, recent determinations of the 
Sun's proper angular momentum Lq = (190.0 it 1.5) x 10^^ kg m^ s~^ from 
helioseismology (Pijpers 1998; 2003), accurate to 0.8%, yield a precessional 
effect one order of magnitude smaller for Mercury (Ciufolini and Wheeler 
1995; lorio 2005c). See Table [T] for the gravitomagnetic precessions of the 
four inner planets. As can be seen, they are of the order of 10~^ — 10^^ " 
cy-^. 

So far, the LT effect on the orbits of the Sun's planets was believed to 
be too small to be detected (Soffel 1989). lorio (2005c) preliminarily in- 
vestigated the possibility of measuring such tiny effects in view of recent 
important developments in the planetary ephemerides generation. It is re- 
markable to note that the currently available estimate of Lq is accurate 
enough to allow, in principle, a genuine test of GTR. Moreover, it was de- 
termined in a relativity-free fashion from astrophysical techniques which do 
not rely on the dynamics of planets in the gravitational field of the Sun. 
Thus, there is no any a priori 'memory' effect of GTR itself in the adopted 
value of Lq . 
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3 Compatibility of the determined extra-precessions 
of planetary perihelia with the LT effect 



3.1 The Keplerian orbital elements 

The Keplerian orbital elements like w are not directly observable quanti- 
ties like right ascensions, declinations, ranges and range-rates which can be 
measured from optical observations, radiometric measurements, meridian 
transits, etc. They can only be computed from a state vector in rectangu- 
lar Cartesian coordinates which also allows to compute predicted values of 
the observations. In this sense, speaking of an "observed" time series of a 
certain Keplerian element would mean that it has been computed from the 
machinery of the data reduction of the real observation^. Keeping this in 
mind, it would be possible, in principle, to extract the LT signal from the 
planetary motions by taking the difference between two suitably computed 
time-series of the Keplerian elements in such a way that it fully accounts 
for the gravitomagnetic signature. Such ephemerides, which should share 
the same initial conditions, would differ in the fact that one would be based 
on the processing of the real data, which are presumed to fully contain also 
the LT signal, and the other one would, instead, be the result of a purely 
numerical propagation. The dynamical force models with which the data 
are to be processed and the numerical ephemeris propagated do not contain 
the gravitomagnetic force itself: only the general relativistic gravitoelectric 
terms must be present. Moreover, the astronomical parameters entering the 
perturbations which can mimic the LT signature should not be fitted in the 
data reduction process: they should be kept fixed to some reference values, 
preferably obtained in a relativity-independent way so to avoid 'imprint- 
ing' effects. Thus, in the resulting "residual" time series AtUobsC*)) the LT 
signature should be entirely present. 

3.2 The EPM2004 ephemerides 

A somewhat analogous procedure was recently implemented with the Ephemerides 
of Planets and the Moon EPM2004 (Pitjeva 2005a; 2005b) produced by 
the Institute of Applied Astronomy (lAA) of the Russian Academy of Sci- 
ences (RAS). They are based on a data set of more than 317,000 observa- 
tions (1913-2003) including radiometric measurements of planets and space- 
fill practice, this procedure cannot be always performed at all times because the obser- 
vations collected at a given time not always allow for the computation of the entire state 
vector in rectangular Cartesian coordinates. 
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craft, astrometric CCD observations of the outer planets and their satehites, 
and meridian and photographic observations. Such ephemerides were con- 
structed by the simultaneous numerical integration of the equations of mo- 
tion for all planets, the Sun, the Moon, 301 largest asteroids, rotations of the 
Earth and the Moon, including the perturbations from the solar quadrupolar 
mass moment J2 and asteroid ring that lies in the ecliptic plane and con- 
sists of the remaining smaller asteroids. In regard to the post-Newtonian 
dynamics, only the gravitoelectric terms, in the harmonic gauge, were in- 
cluded (Newhah et al. 1983). 

3.3 The measured extra-precessions of the planetary perihe- 
Ha and the Lense-Thirring effect 

As a preliminary outlook on the measurability of the Lense-Thirring per- 
ihelion precessions, let us make the following considerations. The magni- 
tude of the gravitomagnetic shift of the Mercury's perihelion over a 90-years 
time span like that covered by the EPM2004 data amounts to 0.0018 ". 
The accuracy in determining the secular motion of Mercury's perihelion can 
be inferred from the results for the components of the eccentricity vector 
k = ecosw and h = esintu reported in Table 4 by Pitjeva (2005b). Indeed, 
the formal standard deviations of k and h are 0.123 and 0.099 milliarcsec- 
onds, respectively. Thus, the formal error in measuring w is about 0.0007 
". An analogous calculation for the Earth yields an error in of 8 x 10~^ 

The EPM2004 ephemerides were used to determine corrections A-djobs 
to the secular precessions of the longitudes of perihelia of the inner planets 
as fitted parameters of a particular solution. In Table 3 by Pitjeva (2005a), 
part of which is reproduced in Table [21 it is possible to find their values 
obtained by comparing the model observations computed using the con- 
structed ephemerides with actual observations. Note that in determining 
such extra-precessions the PPN parameters (Will 1993) 7 and /3 and the so- 
lar even zonal harmonic coefficient J2 were not fitted; they were held fixed 
to their GTR values, i.e. 7 = /3 = 1, and to = 2 x 10"''. Note also that 
the unit values of (3 and 7 were measured in a variety of approaches which 
are independent of the gravitomagentic force itself. Although the original 
purposqj of the determination of such corrections was not the measurement 

*The goal by Pitjeva (2005a) was to make a test of the quality of the previously obtained 
general solution in which certain values of /3, 7 and J2, were obtained. If the construction 
of the ephemerides was satisfactory, very small "residual" effects due to such parameters 
should have been found. She writes: "At present, as a test, we can determine [...] the 
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Table 2: Determined extra-precessions Adjobs of the longitudes of perihe- 
lia of the inner planets, in " cy~^, by using EPM2004 with /3 = 7 = 1, 
= 2 X 10^^. The gravitomagnetic force was not included in the adopted 
dynamical force models. Data taken from Table 3 of (Pitjeva 2005a). It 
is important to note that the quoted uncertainties are not the mere for- 
mal, statistical errors but are realistic in the sense that they were obtained 
from comparison of many different solutions with different sets of parame- 
ters and observations (Pitjeva, private communication 2005a). The correla- 
tions among such determined planetary perihelia rates are very low with a 
maximum of about 20% between Mercury and the Earth (Pitjeva, private 
communication 2005b). 
Mercury Venus Earth Mars 

-0.0036 ± 0.0050 0.53 ±0.30 -0.0002 ± 0.0004 0.0001 ± 0.0005 



of the LT effect, the results of Table 3 by Pitjeva (2005a) can be used to 
take first steps towards an observational corroboration of the existence of 
the solar gravitomagnetic force. Indeed, the uncertainties in the predicted 
values of the LT precessions induced by the error in Lq (Pijpers 1998; 2003) 
amount to 1 x 10~^ " cy~^ for Mercury, 7 x 10"'^ " cy~^ for the Earth and 
2 X 10"'^ " cy~^ for Mars: they are far smaller than the errors in Table EJ so 
that a genuine comparison with the measured precessions make sense. 

By comparing Table [T] and Table [2] of this paper it turns out that 
the predictions of GTR for the LT effect are compatible with the small 
determined corrections to the secular motions of the planetary perihelia 
foii Mercury (-0.0086" cy-^< -0.0020" cy'^< 0.0014 " cy"^), the Earth 
(-0.0006" cy-i< -0.0001" cy-^< 0.0002 " cy-^) and Mars (-0.0004 " cy-^< 
-3x 10~^ " cy~^< 0.0006 " cy~^). In normalized units ^ (/Ugtr = 1) we have 
^Mercury ^ 18^2.5, /i^bT*" = ^±4 and /i^^''^ = -3.3±16. Figure I] summa- 
rizes the obtained results. The discrepancies between the predicted and the 
determined values are reported in Table [3] of this paper. They are smaller 
than the measurement uncertainties, so that a = X] { ^e^ ) ~ ^-^ 

corrections to the motions of the planetary periheUa, which allows us to judge whether the 
values of /3, 7, and J2 used to construct the ephemerides are valid.". The smallness of the 
extra-perihelion precessions found in her particular test-solution is interpreted by Pitjeva 
as follows; "Table 3 shows that the parameters f3 = 1, 7 = 1, and J2 = 2 x 10~^ used to 
construct the EPM2004 ephemerides are in excellent agreement with the observations." 

^In the case of Venus the discrepancy between the predicted and the measured values 
is slightly larger than the measurement error. For such a planet the perihelion is not a 
good observable because of the small eccentricity of its orbit (evenus = 0.0066). 
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Figure 1: The horizontal dash-dotted hnes represent the predicted values 
of the LT secular precessions of the perihelia of Mercury, Venus, the Earth 
and Mars according to GTR. The vertical solid lines represent the values of 
the additional secular precessions of Mercury, Venus, the Earth and Mars 
determined by Pitjeva (2005a) along with their error bars. The predictions 
of the LT effect by GTR (1 in normalized units) are compatible with them 
for Mercury (1.8 ± 2.5), the Earth (2 ± 4) and Mars (-3.3 ± 16.6). 
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Table 3: Comparison between the predicted values (P) of the LT preces- 
sions of the perihelia of Mercury, the Earth and Mars (Table [1]) and the 
determined values (D) of the extra-precessions of their perihelia (Table [2]). 
Their differences are smaller than the errors (E). From them a = 0.2 can 
be obtained. Note that the for zero is 0.8. 

Planet P-D (" cy"^) E (" cy"^) 

Mercury 0.0016 0.0050 

Earth 0.0001 0.0004 

Mars -0.0001 0.0005 



obtained. It must be noted that the determined extra-precessions of Table 
[5] are also compatible with zero, but at a worse level. Indeed, = 0-8 in 
this case. 

A way to improve the robustness and reliability of such a test would 
be to vary the adopted values for the solar oblateness within the currently 
accepted ranges and investigate the changes in the fitted values of the extra- 
precessions. Moreover, it would also be important to produce an analogous 
set of solutions with /?, 7 and Jg* fixed in which the extra-precessions of the 
nodes are determined; in this way it would be possible to use only Mercury. 



3.4 Some possible systematic errors due to other competing 
effects 

In order to check our conclusion that the LT effect is the main responsible 
for the observed secular corrections to the planetary perihelia Atzjobs let us 
focus on Mercury and on the known perturbations which could induce a 
secular extra-perihelion advance due to their mismodelling. 

The major sources of secular advances of the perihelia are the Schwarzschild 
gravitoelectric part of the solar gravitational field and the quadrupolar mass 
moment J2 of the Sun. Their nominal effects on the longitudes of perihelion 
of the inner planets are quoted in Table H] and Table [5] of this paper; the 
analytical expressions are 

3 nJ2 f ^ 



where n = ^ GMja? is the Keplerian mean motion and R is the mean 
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Table 4: Nominal values of the secular post-Newtonian gravitoelectric pre- 
cessions of the longitudes of perihelion w of Mercury, Venus, Earth and Mars 
in " cy"-*^. Their mismodelled amplitudes are fixed by the uncertainties in 
7 and /3 which are of the order of 0.01% (Pitjeva 2005a). However, devia- 
tions from their relativistic values are expected only at a 10~® — 10~^ level 
(Damour and Nordtvedt 1993). 

Mercury Venus Earth Mars 

42.9812 8.6247 3.8387 1.3509 



Tabic 5: Nominal values of the classical secular precessions of the longitudes 
of perihelion vc of Mercury, Venus, Earth and Mars, in " cy~^, induced by 
the solar quadrupolar mass moment Jf'. The value J2 = 2 x 10~^ used in 
(Pitjeva 2005a) has been adopted. Their mismodelled amplitudes are fixed 
by the uncertainty in J® which is of the order of ~ 10%. 

Mercury Venus Earth Mars 

0.0254 0.0026 0.0008 0.0002 



equatorial radius of the central body. In view of their large size with 
respect to the LT effect, one could legitimately ask if the determined extra- 
precessions arc due to the systematic errors in such competing secular rates. 
An a-priori analytical analysis shows that it should not be the case. 

3.4.1 The impact of the solar oblateness 

In regard to J^^, which is still rather poorly known, only values measured in 
such a way that no a priori 'imprinting' effects occurred should be considered 
for our purposes. 

E.g., the most recent determinations of the solar oblateness based on 
astrophysical techniques yield values close to 2.2 x 10"'' (Paterno et al. 
1996; Pijpers 1998; Mecheri et al. 2004) with discrepancies between the 
various best estimates of the same order of magnitude of their errors, i.e. ~ 
10"^. Let us see if such determinations are compatible with the determined 
extra-advances of perihelia. By assuming a correction of ~ 10% of the 
adopted reference value by Pitjeva, i,e. Jjf^ = 2 x 10~^, the resulting 
residual precession due to the solar oblateness would amount to -1-0.0025" 
cy~^ for Mercury. It falls outside the measured range. 

A way to a priori cancel out any possible impact of the uncertainty in 
the solar oblateness consists in suitably combining the perihelia advances 
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of two planets so to de-correlate by construction the LT effect and the pre- 
cessions due to J2. This approach allows to extract the gravitomagnetic 
signal independently of the solar quadrupolar mass moment. On the other 
hand, it is also possible to measure a correction to it independently of 
the LT effect. In turn, such value for 6J2 can be used to check if it can 
accommodate the determined extra-precessions of Table [2j Let us write 

f A ■ Mercury • Mercury r t0 i • Mercury 

< 

AT^Earth _ ■ Earthr t0 , -Earth 

By solving eq. ([6]) with respect to /u it is possible to obtain 



(6) 



A - Mercury , „ a - Earth 
^ ^^^ercury ^ ^^^Earth ' ^ > 



with 



■ Mercury / Earth \ 7/2 / -, _ 2 \ 

ci = -^W--(i=;) =-30.1930, (8) 



—Earth I ^Mercury / I 1 o2 

^.2 V« V yJ^-^Mercury 



and 



■ Mercury -Earth n nm Q "^,,-1 ln\ 

■ojj^rp + ciWuj^ = O.OOio cy . (y) 

The combination of eq. ([7]) is not affected by the solar oblateness whatever 
its real value is: indeed, 

^Mercury ^ ^^^Earth ^ y (^g) 

By inserting the value of eq. ([5|) , the figures of Table [T] and the results of 
Table [2] in eq. 1^ one obtains for such a combinatiorH ^ohs = 1-8 ± 10. As 
can be noted, the best estimate for /i does not change with respect to the 
case of Mercury's perihelion only, as if departures of the solar oblateness 
from the adopted reference value were of little importance. Indeed, if we 
solve eq. ^ with respect to the correction to the Sun's quadrupolar mass 
moment the equation 

Mercury , ArhEarth 
^jQ ^ ^^obs +^l^^°bs (11) 



®The error has been evaluated by summing in quadrature the errors of Table[2]according 
to eq. 0. 
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with 



■ Mercury / \ 3 / i _ 2 \ 3/2 

di = -^i^— / ^Earth =_18.3864, (12) 



is obtained. Eq. (fTTI) allows to measure the correction to the adopted value 
of J2', by construction, independently of the LT effect in the sense that 

^Mercury ^ ^^^Earth ^ y (^3^ 

The result ifl 

SJ® = (+0.01 ± 0.47) X 10"^. (14) 

Such value can be considered as a dynamical measurement of the solar 
oblateness independent of the general relativistic gravitomagnetic features 
of motion. It induces a "residual" precession of +0.0001 " cy~^ on Mer- 
cury's perihelion, which is smaller than its observed extra-advance and the 
related error. For the Earth the "residual" effect of eq. (I14p would amount 
to +4 X 10"*^ " cy-i. 

It is important to note that the combination of eq. d?]) yields = 0.007, 
while, by assuming zero extra-precessions, one has = 0.03: also in this 
case, the relativistic prediction of LT is in better agreement with data than 
the zero-effect hypothesis. 



3.4.2 The post-Newtonian gravitoelectric precessions 

Although the large nominal values of their precessions, the post-Newtonian 
gravitoelectric terms do not represent a problem. Indeed, they are fully 
included in the dynamical force models of EPM2004 in terms of the PPN 
parameters P and 7 which are presently known at a 10~^ — 10~^ level (Pitjeva 
2005a; Bertotti et al. 2003). Moreover, theoretical deviations from the GTR 
values are expected at a 10~^ — 10~^ level (Damour and Nordtvedt 1993). 



3.4.3 The impact of the asteroids 

As already noted, the dynamical force models adopted in EPM2004 also 
include the action of the major asteroids and of the ecliptic ring which 
accounts for the other minor bodies. Indeed, it has recently pointed out 

^Note that Pitjeva (2005a) derived a negative correction of order 10~* from the 
determined extra-advance of Mercury's perihehon only, without taking into account the 
biasing impact of the LT effect which amounts to ~ 8% for Mercury, as can be inferred 
from Table [Hand Table El 
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that their impact hmits the accuracy of the inner planets' ephemerides over 
time-scales of a few decades (Standish and Fienga 2002) in view of the 
relatively high uncertainty in their masses (Krasinsky et al. 2002; Pitjeva 
2005b) Recently, Fienga and Simon (2005) have shown that also Mercury's 
orbit is affected to a detectable level by secular perturbations due to the 
most important asteroids. 

May it happen that the mismodelled part of such secular precessions 
could explain the observed Aw^^^^^^^l 

From Table 3 by Fienga and Simon (2005) the nominal amplitude of 
the secular perturbations on ^Mercury 295 major asteroids can be 

calculated. It turns out to be 0.0004 " cy~^; even assuming a conservative 
~ 10% uncertainty (Pitjeva 2005b), it is clear that the asteroids are not the 
cause of the determined extra-shift of Mercury's perihelion. 



3.4.4 The impact of non-Einsteinian effects 

In regard to other possible sources of extra-secular precessions of the plan- 
etary perihelia outside the scheme of the Newton-Einstein gravity, recently 
it has been shown by Lue and Starkman (2003) that the multidimensional 
braneworld gravity model by Dvali, Gabadadze and Porrati (2000) predicts 
also a secular perihelion shift in addition to certain cosmological features. 
By postulating that the current cosmic acceleration is entirely caused by the 
late-time self-acceleration, constraints from Type lA Supernovae data yield 
a value of ~ 0.0005 " cy~^ for the Lue-Starkman planetary precessions. Also 
this effect is too small to accommodate the determined additional perihelion 
advance of Mercury. 



3.5 Analysis of other independent data 

The shift of the perihelion yields a variation of the planet's range that can 
be expressed as (Nordtvedt 2000) Ar = eaAzu. In the case of Mercury the 
centennial variation due to the Lense-Thirring precession amounts to -115 
m. 

Fienga et al. (2005) used their numerical ephemeride^ INPOP, recently 
produced at the Institute of mecanique celeste et de calcul des ephemerides 
(IMCCE), to fit different kinds of observations including also the radar- 
ranging data to the inner planets. They found for the Mercury's range 

*They include a complete suit of improved dynamical models, apart from just the post- 
Newtonian gravitomagnetic forces whose effects are, thus, fully present in the determined 
residuals. 
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residual the value Jr, 



meas 



95.6 lb 784 m. It is worth noting that in 



obtaining these results also the impact of the asteroids on Mercury's orbital 
motion was accounted for. As it can be noted, also in this case the errors 
are large, but the general relativistic prediction for the Lense-Thirring effect 
are in better agreement with the data than the hypothesis of null effect 
([(P - D)/E]2 = 6 X 10"^ and [(P - D)/E]2 = 1 x lO^^, respectively). 

4 Constraints on a Yukawa-like fifth force 

The differences between the determined extra-precessions and the predicted 
LT rates of Table [3] of this paper can also be used to strongly constrain, at 
planetary length-scales 10^'' — 10^^ m, departures from the inverse-square- 
law phenomenologically parameterized in terms of the magnitude |a| of the 
strength of a Yukawa-like fifth force (Fischbach et al. 1986; Adelberger et 
al. 2003). Indeed, a potential 



where A is the range of such a hypothesized force, can produce a secular 
perihelion advance over scales A comparable to a (Lucchesi 2003) 



By using the figures of Table [3] it is possible to constrain aio"^ 10^^^ — 10^^^ 
level atr = A~lA.U. The most recently published constraints in the 
planetary range are at 10^^ — 10"^'' level (Bertolami and Paramos 2005; 
Reynaud and Jaekel 2005). 

5 Discussion and conclusions 

In this paper we discussed the possibility of performing new tests of post- 
Newtonian gravity in the Solar System. To this aim, we analyzed the esti- 
mated corrections to the secular rates of the perihelia of the inner planets 
of the Solar System recently determined by E.V. Pitjeva (Institute of Ap- 
plied Astronomy, Russian Academy of Sciences). She used the EPM2004 
ephemerides with a wide range of observational data spanning almost one 
century; in a particular solution, she solved also for the secular motions of 
the perihelia by keeping fixed the PPN parameters (3 and 7 and the solar 
quadrupole mass moment J^, and neglecting the gravitomagnetic force in 
the dynamical force models. 




(15) 




(16) 
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It turns out that the post-Newtonian LT secular precessions predicted 

by GTR are compatible with the determined extra-prcccssions for Mercury, 
the Earth and, to a lesser extent. Mars: in normalized units (/i = 1 in GTR) 
we have ii^hs = 1-8 ± 2.5 for Mercury, //obg = 2 it 4 for the Earth and n^hs = 
—3.3 =b 16.6 for Mars. A suitable combination of the perihelia of Mercury 
and the Earth, which cancels out any possible bias by J^, yields /Uobs = 
1.8 ± 10. It must be noted that the errors are still large and the data are 
compatible also with the hypothesis of zero extra-precessions, but at a worse 
level with respect to the relativistic LT prediction. If confirmed by further, 
more extensive and robust data analysis by determining, e.g., the extra- 
precessions of the nodes as well, it would be the first observational evidence 
of the solar gravitomagnetic field. The processing of further amounts of 
data, along with those expected in future from the forthcoming planetary 
mission BepiColombo and, perhaps. Messenger and Venus Express as well, 
although to a lesser extent, will further improve the accuracy in determining 
the orbital motion of these planets and, consequently, the precision of the 
LT tests. 

A by-product of the present analysis is represented by new, strong con- 
straints (10~^^ — 10~^^) on the strength of a Yukawa-like fifth force at scales 
of about one Astronomical Unit. 
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